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Abstract 

On the basis of both a conventional relativistic nuclear fluid dynaaical aodel and a 
two- fluid generalization that takes into account the interpenetration of the target and 
projectile upon contact, we calculate collisions between heavy nuclei aoving at relativis- 
tic speeds. This is done by solving the relevant equations of aotion nuaerically in three 
spatial diaensions by use of particle-in-cell finite-differe;.^e coaputing techniques. We 
study the effect of incorporating a density isoaer, or quasistable state, in the nuclear 
equation of state ac three tiaes noraal nuclear density, as f^ll as the effect of doubling 
the nuclear coapressibility coefficient. For the reaction **Ne + ***U at a laboratory boa 
bar.:ing energy per nucleon of 393 HeV, the calculated distributions in energy and angle of 
outgoing charged particles are compared with recent experimental data both integrated over 
all impact parameters and for nearly central collisions. 

Introduc tion 

Our colloqui«im on drops and bubbles has brournt togdthev people in such diverse disci- 
plines as astrophysics, fluid dynamics, and nuclear physics, working with physical systems 
ranging in size from stars through water drops to atomic nuclei. Having already heard 
from astrophysicists and fluid dynamicists, it is time now to turn our attention from some 
of the largest objects in the universe to some of the smallest. 

An atomic nucleus, with diameter about 10*'* m and mass about 10'^^ kg, is >.bout 10‘* 
times as dense as water. Discovered in 1911 by Sir Ernest Rutherford, it is to lowest or- 
der a collection of positively charged protons .-.r.-i uncharged neutrons held tightly togeth- 
er by short-range nuclear forces. As first suggested by Niels Bohr in 1936, a nucleus be- 
haves in some respects like a uniformly charged drop of liquid with surface tension. 

The liquid-drop model has been widely used to explain the breakup of a heavy nucleus 
into two smaller fragments in the process of nuclear fission, as well as to explain many 
features associated with collisions between two nuclei. Until recently, the bombarding 
energy in such collisions has been fairly low, with the result that the nuclear density 
remains close to its equilibrium value and the excitation energy is relatively low. How- 
ever, a few years ago accelerators were developed that can accelerate heavy nuclei to rel- 
ativistic speeds, and we are now beginning to explore what happens when nuclei become high 
ly compressed and excited. 

The collision of heavy nuclei at high energy is extremely complicated, lying somewhere 
in between two limiting possibilities that have been used to describe it. At one extreme, 
the process could be dominated by a series of collisions between individual particles mak- 
ing up the nuclei or produced in the reaction. This first limit, which would be realized 
if the particle mean free path were much longer than the nuclear force range, has been 
studied in terms of microscopic approaches such as the relativistic intranuclear cascade, 
where the basic input is experimentally measured two-particle cross sections. 

At the other extreme, the process could be dominated by coherent collective-field ef- 
fects, resembling instead the collision of two drops of nuclear fluid. This second limit 
would be realized if there are many degrees of freedom, sufficient time during the colli- 
sion to establish local equili'onum, and a short mean free path for stopping a particle. 

In relativistic nuclear collisions of the type considered here, the first two conditions 
are satisfied moderately well whereas the last condition is more uncertain. Estimates 
based on collisions between two particles in free space give relatively long mean free 
paths, but the mean free path could be reduced significantly inside a nuclear medium be- 
cause of many-body effects or a nuclear phase transition.* This second limit has been 
studied by several different groups around the world in terms of fluid dynamical models, 
where the basic input is the nuclear equation of state. Ve will be concentrating here on 
recent work perform -.d at Los Alamos within this approach. 



Nuclear equation of state 

One reason for studying relativistic nuclear collisions is to learn about the nuclear 
equation of state, the fundaaental relationship specifying how pressure depends upon den- 
sity and themal energy. At present we have experiaental information about this important 
function only in the vicinity of the equilibrium ground state. However, theoretical spec- 
ulations suggest that it may be extremely complicated, with nuclear matter undergoing one 
or more phase transitions as its density is increased.^ This is illustrated schematically 
in Figure 1, which shows how the ground-state energy per nucleon Eo(n), or zero- tempera- 
ture compressional energy per nucleon, might depend upon nucleon number density n. We 
know experimentally that the equilibrium energy per nucleon Eo(i>o) = *1^ HeV, that the 
equilibrium density n,> = 1.5 x 10** nucleons/ar , and that the nuclear compressibility co- 
efficient K = 210 MeV. What happens away from equilibrium is currently unknown, but dou- 
bling the nuclear density from its normal value could lead to a pion condensate, or a 
state containing a large number of bound pions. Compression to several times normal densi- 
ty could result in a density iosmer, or a quasistable state existing at other than normal 
density. Still further compression could produce quark matter, in which the quarks that 
comprise nucleons become free. To determine whether or not any of these phase transitions 
actually exist in nuclei is the exciting challenge that we face! 


In addition to the above compressional energy, nuclear matter at rest can contain ther- 
mal energy, so that the total internal energy per nucleon is 

E(n,I) = fo(n) + I , (1) 


where I is the thermal energy per nucleon, which is itself a function of n and either the 
entropy per nucleon S or the temperature T. The pressure p is then given by*’’ 
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containing separate contributions from the compressional energy and the thermal energy. 
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Figure 1. Compressional energy, illus- 
trating three conjectured nuclear phase 
transitions . 



Nucleon Number Density n/n^ 

Figure 2. Compressional energy in three 
of our nuclear equations of state. 
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Foi the grcund-state energy per nucleon E<>(n) we use both an old functional fora^>* 

Eo(n> = a n^'^^ - b n + c n^'^^ , (3) 

with positive constants a. b, and c that yield a compressibility coefficient K = 294.8 
HeV, and a new functional form^ producing the three curves sho%m in Figure 2. The solid 
line shows the result for compressibility coefficient K = 200 MeV, and the dashed line 
shows the effect of doubling the compressibility coefficient to 400 MeV. The dot-dashed 
curve shows the result for a density isomer at a density that is three times normal nu- 
clear density, with an energy 2 HeV higher than at normal density and with the same cur- 
vature. In each of these three cases we use the value Eo(no) - -8 MeV to simulate the 
loss in binding energy for finite nuclei arising from surface and Coulomb energies. The 
new func tional form has the property that the speed of sound approaches the speed of light 
in the imit of infinite compression. This is achieved by parametrizing Eo(n) for n great- 
er than a critical value in terms of three smoothly joined parabolas in the square root 
of the density, so that in the limit of infinite compression it increases linearly with 
density. 

For the thermal contribution to the pressure we use the result*'^ 

Pthe™! = Isis = 3 ^ 

which is derived from the nonrelativistic Fenni-gas model. 

Relativistic nuclear fluid dynamics 

In a complete nuclear fluid dynamical calculation, we would need to take into account 
nuclear energy. Coulomb energy, nuclear viscosity, thermal conductivity, and single-par- 
ticle effects, as well a-^ the production of additional particles and the associated radia- 
tive loss of energy from the system. However, in nuclear collisions of the type consider- 
ed here, these effects are small compared to those caused by the dominant kinetic, compres- 
sional, and thermal energies, and are consequently neglected. 

The covariant relativistic fluid dynamical equations that we solve express the conserva- 
tion of nucleon number, momentum, and energy, for a specified nuclear equation of state. 

In units in which the speed of light c = 1, these equations are®‘^ 


II ♦ V-(On) = 0 , (5) 

II + y. (Ofi) = -Vp , (6) 

and 

II + V-(Oe) = -V-(Op) . (7) 

where N, ft, and E are respectively the nucleon number density, momentum density, and ener- 
gy density (including rest energy) in the laboratory reference frame and v is the velocity 
of matter relative to the laboratory frame. The three laboratory- frame quantities are re- 
lated to rest-frame quantities by 


N = yn, (8) 

ft = Y^(e + p)v , (9) 

and 

E = y^(e + p) - p , (10) 

where y = (1 - v*)'^ and e is the internal energy density in the rest frame, which is re- 
lated to the internal energy per nucleon of Eq. (1) by 


e = 1 mo + E( n , I ) 1 n 


( 11 ) 


including the nucleon mass mo. 
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For a given nuclear equation of state and for given initial conditions, we solve these 
equations nuaerically in three spatial diaeasions by use of a particle-in-cell finite-dif- 
ference coaputing aethod.^'* In this technique, the fluid is represented by discrete La- 
grangian coaputational particles, which aove through a aesh consisting of fixed cubical 
Eulerian cells. Froa finite-difference representations of Eqs. (5)-(7), the values of N, 

H, and E for (ach cell are calculated at later tiaes in teras of preceding values. The 
values of n, v, e, and p throughout the aesh are obtained by aeans of a partial algebraic 
reduction followed by the iterative solution of a transcendental equation in one unknown. 

Soae exaaples of the solutions are shown in Figure 3 for the reaction *°Ne *■ at a 

laboratory b^barding energy per nucleon of 250 HeV, corresponding to an incident speed 
that is S2X the speed of light. The nuclear equation of state is given by Eqs. (2)-(4). 
Each coluan presents a side view of the aatter distribution evolving in tiae for a differ- 
ent iapact paraaeter. The initial fraae in each case shows a target boabarded froa 

above by a Lorentz-contracted *®Ne projectile. The projectile and target are represented 
by coaputational particles, which are initially aligned so that in the direction perpendic- 
ular to the page only a single point is visible. However, as the iapulse resulting froa 
the collison propagates throughout the systea this alignaent is destroyed and additional 
particles coae into view. 

The characteristic features of the tiae evolution vary sysceaatically with iapact paraa- 
eter. In nearly central collisions the target and projectile are substantially deforaed, 
coapressed, and excited, with curved shock waves produced. These are followed by rarefac- 
tion waves and an overall expansion of the natter into a aoderately wide distribution of 
angles. At the other extreae, in peripheral collisions the projectile is fragaented into 
a portion that proceeds roughly straight ahead at its original velocity and another por- 
tion that deposits its energy in the target. This disturbs the target auch less violently 
than in nearly central collisions, and its deformation, compression, and excitation are 
therefore much less. 


We show in Figure 4 the effect of varying the nuclear equation of state on nearly cen- 
tral collisions at the higher laboratory bombarding energy per nucleon of 393 HeV, corre- 
sponding to an incident speed that is 71% the speed of light. The results for the differ- 
ent equations of state are very similar to one another, but for our equation of state with 
a density isoaer the expansion starts somewhat later because the natter is coapressed to a 
higher density than for our two conventional equations of state. 
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Figure 3. Tine evolution of the natter 
distribution for three impact parameters. 
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Figure 4. Tine evolution of the natter 
distribution for three nuclear equations 
of state. 
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For our conventional equation of state with compressibility coefficient K = 200 HeV, 
the matter is compressed to a maximum rest-frame density of 3.9 no and remains above 2 n<> 
for 1.0 X 10'^^ s. For our stiffer equation of state with K = 400 HeV the matter is com- 
pressed to a meucimum rest-frame density of 3.4 n^ and remains above 2 n<> for 0.7 x 10'** s. 
For our softer equation of state with a density isomer the matter is compressed to a maxi- 
mum rest-frame density of 5.1 no and remains above 3 no for 1.4 x 10'** s and above 2 no 
for 1.9 X 10'** s. 


Comparison with experimental data 

For the reaction *®Ne ♦ ***U at a laboratory bombarding energy per nucleon of 393 HeV, 
Sandoval et al.* have recently measured at the Berkeley Bevalac the distributions in ener- 
gy and angle of outgoing charged particles, including contributions from protons, deuter- 
ons, tritons, *He particles, and ♦He particles. Also, by selecting only those events con- 
taining large numbers of associated charged particles. Stock et al.*^ have measured for 
nearly central collisions the distributions in energy and angle of outgoing charged par- 
ticles, including contributions from protons, deuterons, and tritons. Unfortunately, this 
selection does not eliminate possible contributions from large impact parameters, for 
which the high-energy yield in forward directions is larger than for central collisions. 
These experimental data in the form of energy spectra at four laboratory angles ranging 
from 30° to 150° are shown by the solid circles in Figure 5, where they are compared with 
histograms calculated for our conventional equation of state with compressibility coeffi- 
cient K = 200 NeV. We have also made similar comparisons with results calculated for the 
other two equations of state illustrated in Figure 2. 

Our calculated distributions are obtained by constructing from the velocity vectors at 
some large time the energy and angular distributions for the expanding matter. The small 
amount of matter that already has passed through the top and side boundaries of the com- 
putational mesh is also included. By integrating over the appropriate ranges of impact 
parameter, we compute the double-differential cross section corresponding both to all im- 
pact parameters and to central collisions constituting 15% of the total cross section. 

The cross section for the outgoing matter distribution is then converted into the cross 
section d*o/dEdU for outgoing charged particles under the assumption of uniform charge den- 
sity. Some measure of the fairly large numerical inaccuracies inherent in fluid dynamical 
calculations can be determined from the fluctuations in the histograms, which are obtained 
using angular bins of 10° width. 

Examining first the results for all impact parameters given in left-hand side of Figure 
5, we see that at low energy the calculated results are for all angles higher than the ex- 
perimental results. This is because of our neglect of binding, which causes the entire 
system to completely disintegrate into slowly moving matter for an arbitrarily small im- 
pulse. A* higher energy the calculations with all equations of state reproduce, to within 
numeric uncertainties, the experimental data at all angles. We have found that the nu- 
clear equation of state has little effect on the single-pa-ticle-inclusive cross section 
d*o/dEdU integrated over all impact parameters. 

We turn now to the results for central collisions, which are given in the right-hand 
side ot Figure 5 for our conventional equation of state with K = 200 MeV. At intermediate 
angles the results calculated with the three equations of state are very similar to one 
another, to within numerical uncertainties. However, at 6 = 30° the slope of the energy 
spectrum decreases significantly as we go from a stiff equation of state with K = 400 MeV 
through an intermediate one with K = 200 HeV to a soft one that contains a density isomer. 
Also, at t - 150° the results calculated for the density isomer are somewhat higher than 
those calc. lated for the two conventional equations of state. These differences arise be- 
c.iuse rhe softer density-isomer equation of state leads to higher initial density and ther- 
mal ent igy per nucleon, which increases the thermal contribution to the cross section in 
reg ...IS where it would otherwise be small. 

Because of our neglect of binding, at low energy the calculated results for central col- 
lisions are also higher than the experimental results for all angles and equations of 
state except for 8 = 150° with the two conventional equations of state, where the rapid 
expansion in the backward direction suppresses the cross section. At higher energy the 
calculations with all equations of state reproduce, to within numerical uncertainties, the 
i!xperimental data at all angles except 6 = 30°, where the calculated energy spectra for 
oill three equations of state have significantly larger slopes than the experimental spec- 
trum. Ihis important discrepancy for central collisions in the forward direction could 
arise from several different possibilities, but the most likely is that upon contact the 
ra' <;et and projectile interpenetrate substantially. This interpenetration can be taken 
' nto account, while retaining some degree of coherent collective flow, by means of a nu- 
clear two-fluid model, to which we now turn our attention. 
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Nuclear two- fluid dynaaics 

In Che nuclear two- fluid nodel coupled relativistic equations of notion are solved for 
separate target and projectile nuclear fluids.**** The terns in the equations that couple 
the two nuclear fluids are obtained fron the cross section and nean longitudinal nonentun 
transfer for free nucleon-nucleon collisions. At low relative velocities the target and 
projectile fluids nerge, in which case conventional relativistic nuclear fluid dynanics 
(which we alternatively refer to as one-fluid dynanics) is recovered. 


The equations of notion for the target and projectile nuclear fluids express the con- 
servation of nucleon nunber, energy, and nonentun, plus the transfer of energy and nonen- 
tun between me two fluids. For a given fluid, the relativistic equations of notion are 
analogous to Eqs. (5)-(7) for conventional nuclear fluid dynanics, but contain in addition 
coupling terns that describe the transfer of energy and nonentun fron one fluid to the 
other as they interpenetrate. In particular, when the relative velocity of the two fluids 
is large conpared to the Femi velocity of the nucleons conprising each fluid, the rela- 
tivistic equations of notion for fluid 1 are*’** 


♦ V.(VjNj) = 0 , (12) 


3M, 

♦ V'(vjftj) = -Vp^ - D(yjVj - > 2 ^ 2 ) * (13) 

and 

3E, 

♦ V-(vjEj) = -V-(VjPj) - D(ifj - Yj) ; (14) 

those for fluid 2 are obtained by interchanging the subscripts 1 and 2. The drag function 
D involves the nucleon nunber densities and velocities of the two fluids, as well as the 
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Figure 5. Conparison of calculated histo- 
grans with experinental points (Refs. 9 
and 10). 
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Figure 6. Tine evolution of the natter 
distriubtion calculated in the two-fluid 
node I . 



cross section and average longitudinal aomentum transfer for free nucleon-nucleon colli- 
sions. For the nuclear equation of state specified by Eqs. (2)-(4), these equation of mo 
tion are solved by use of a two-fluid generalization of the relativistic particle-in-cell 
technique.®* •> ** 

We show in Figure 6 some examples of the two-fluid solutions for the reaction *®Ne + 
at a laboratory bombarding energy per nucleon of 250 MeV. These results are quali- 
tatively similar to those shown in Figure 3 for conventional nuclear fluid dynamics, es- 
pecially during the later stages of the process. However, because the target and projec- 
tile interpenetrate somewhat before they begin to respond to the presence of the other, 
the initial compression is less and the amount of matter emerging in the backward direc- 
tion is less in the two- fluid model than in conventional nuclear fluid dynamics. 


Figure 7 shows the resulting energy spectra calculated for central collisions in both 
the one- fluid and two- fluid models with the same equation of state, along with the experi- 
mental data of Stock et al.^® For the three angles 6 = 70“ , 110“, and 150“, the two mod- 
els reproduce equally well, to within numerical uncertainties, the experimental data at 
high energy and are both larger than the experimental data at low .energy because of our 
neglect of binding. For 6 = 30“ ; the two-fluid model agrees with the experimental data 
substantially better than does the one-fluid model, although the slope calculated from the 
two- fluid model is still somewhat, larger than the experimental slope. However, as shown 
earlier for the one-fluid model, a softer nuclear equation of state than the one used here 
would increase slightly the calculated high-energy yield in forward directions. 


An alternative and perhaps more illuminating way of making the comparisons for central 
collisions is in the form of angular distributions for fixed outgoing laboratory momentum 
per nucleon, as shown in Figure 8. The experimental angular distributions for low out- 
going momenta contain a small peak that shifts to smaller angles and finally disappears 
for higher outgoing momenta Unfortunately, this peak could be the result of either the 
neglect of ®He, ^He, and heavier composite particles, whose yields are concentrated at low 
energies and forward directions, or the Coulomb interaction, which provides a transverse 
driving force for slow-moving charged particles. 
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Figure 7. Comparison of calculated histo- Figure 8. Comparison of calculated histo 

grams with experimental points (Ref. 10). grams with experimental points (Ref. 10). 
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The histograms are calculated with energy bins of 20 MeV. Because of our neglect of 
binding, the results for the lowest outgoing momentum calculated in both the one-fluid and 
two-fluid models are much larger than the experimental data. With increasing outgoing mo- 
mentum, the one- fluid model predicts angular distributions that are narrower than the ex- 
perimental distributions and that are peaked at increasingly larger angles, which is op- 
posite to the experimental trend. However, as shown by St&cker et al.,^ the agreement be- 
tween the one- fluid model and experiment would be improved by superimposing the remaining 
thermal energy at a freezeout density at which fluid dynamics ceases to be valid. The two- 
fluid model predicts angular distributions with peaks that shift to smaller angles with 
increasing outgoing momentum, as is observed experimentally. For intermediate outgoing 
momenta the experimental data are lower in absolute value than the two- fluid calculations, 
but for the two highest outgoing momenta the experimental data agree with the two- fluid 
calculations to within their numerical uncertainities. By comparison, microscopic models 
based on independent two-particle collisions yield for central collisions angular distribu- 
tions that are essentially forward peaked at all outgoing energy. 

Conclusions 


With some qualifications, we conclude that in relativistic nuclear collisions the tar- 
get and projectile interpenetrate substantially upon contact, but that some degree of co- 
herent collective flow is involved. Although conventional relativistic nuclear fluid dy- 
namics is deficient in several respects, relativisitic nuclear two-fluid dynamics satis- 
factorily describes many aspects of relativistic nuclear collisions. Because the calculat- 
ed results are not very sensitive to the input nuclear equation of state, we do not yet 
know whether or not there are any phase transitions as the nuclear density is increased. 

The answer to this question is relevant not only to atomic nuclei but also to neutron stars, 
illustrating once again how drops and bubbles help unify our understanding of physical sys- 
tems differing in size by twenty orders of magnitude. 
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